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We  present  experimental  data  for  BaTi03  that  exhibit  an  intensity  dependence  in  the  absorption  coefficient  and  the 
two-beam  coupling  coefficient  at  intensities  between  0.002  and  40  W/em2'  The  effective  empty -trap  concentration 
was  found  to  increase  with  intensity.  We  present  a  model,  in  the  spirit  of  commonly  used  (photorefractive)  theories, 
for  photorefraction  and  optical  absorption  that  explains  these  effects.  The  intensity  dependence  was  attributed  to 
the  presence  of  secondary  photorefractive  centers,  t ' 


* , 


INTRODUCTION 

Barium  titanate  (BaTiO;))  has  been  the  subject  of  extensive 
research  in  recent  years.  Through  the  photorefractive  ef¬ 
fect,  it  is  able  to  show  high  optical  nonlinearities  at  low 
power  levels  and  has  been  used  for  applications  such  as 
optical  phase  conjugation1  and  optical  resonators.23  The 
physical  process  involved  in  the  photorefractive  effect  is  a 
light-induced  redistribution  of  charges  among  traps.  The 
resulting  space-charge  field  then  modulates  the  refractive 
index  through  the  linear  electro-optic  effect. 

Although  there  has  been  much  interest  in  the  applications 
of  the  photorefractive  effect,  little  is  known  regarding  the 
species  responsible  for  the  effect  in  BaTi03.  Previous  re¬ 
search  suggests  that  iron  impurities  are  the  primary  sources 
and  traps  of  the  photorefractive  charge  carriers.4-5  Models 
based  on  these  deep  traps  provide  a  reasonably  good  under¬ 
standing  of  the  photorefractive  effect  but  are  yet  incom¬ 
plete.6-9  For  example.  Motes  and  Kim  found  an  intensity 
dependence  in  the  absorption  and  two-beam  coupling  coeffi¬ 
cients  that  is  not  explained  by  the  current  models.1011 

In  this  paper  we  examine  the  intensity-dependent  absorp¬ 
tion  and  photorefractive  effects  in  BaTiOa  and  present  a 
model  in  which  the  absorption  characteristics  are  accounted 
for  by  the  presence  of  secondary  centers.  These  are  inter¬ 
mediate-level  traps  that  are  highly  ionized  at  room  tempera¬ 
ture  and  acquire  free-charge  carriers  generated  by  the  inter¬ 
action  of  the  laser  light  with  deep-level  impurities.  We 
incorporate  the  secondary  centers  into  a  photorefractive 
model  and  show  that  they  lead  to  intensity-dependent  pho¬ 
torefractive  effects  not  predicted  by  current  models.  Char¬ 
acteristics  of  these  centers  are  discussed. 

GENERAL  BACKGROUND 

The  grating  formation  in  a  photorefractive  material  is  com¬ 
monly  described  by  either  the  band -transport  model6-8  or 
the  hopping  model,  both  of  which  give  the  same  results. 
The  band  model  is  shown  schematically  in  Fig.  1.  The 
charge  carriers  may  be  either  electrons  or  holes  (or  both). 
Phot...  onductivity  in  BaTiO;i  is  usually  hole  dominated, 
and,  for  simplicity,  we  discuss  the  single-carrier  model.  X 


and  X'  constitute  a  complete  set  of  recombination  centers. 
Holes  are  photoionized  from  X  and  recombine  at  X'.  Over¬ 
all  charge  neutrality  is  maintained  by  other  states  that  are 
not  photorefractively  active.  Under  illumination  by  sinus¬ 
oidally  modulated  light,  as  from  an  interference  pattern 
between  two  coherent  optical  beams,  holes  optically  excited 
in  the  high-intensity  regions  migrate  by  diffusion  and  drift 
to  the  darker  regions  of  the  crystal,  where  they  recombine 
with  empty  traps.  This  redistribution  of  charges  results  in 
an  ionic  space-charge  grating  that  is  out  of  phase  with  the 
periodic  light- intensity  pattern.  The  charge  separation  is 
balanced  by  a  periodic  space-charge  electric  field  that  is  jt/2 
phase  shifted  with  respect  to  the  light  intensity  (no  external¬ 
ly  applied  electric  field).  This  space-charge  field  then  mod¬ 
ulates  the  refractive  index  through  the  electro-optic  effect. 
The  phase  grating  is  also  jt/2  shifted  with  respect  to  the  light 
intensity. 

A  mathematical  description  of  the  grating-formation  pro¬ 
cess  has  been  given  by  Kukhtarev  et  al .6  for  a  light-intensity 
pattern  resulting  from  the  interference  of  two  coherent 
beams  Is  and  Ip: 

Hz)  =  /0[1  +  m  exp(iKz)],  (1) 

where  m  is  the  modulation  index  given  by  m  =  2(lslp)l/2  /«,  /o 
=  Is  +  Ip,  and  K  is  the  grating  wave  number.  For  small 
modulation  index,  the  magnitude  of  the  modulated  space- 
charge  field  can  be  written  as4-5 

£sc  *  R(mkBT/q)\K/[\  +  (X/X0)2] j,  (2) 

where 

K02  =  q2NE/(((okBT),  (3) 

Ne  =  Nxil  -  Nx/Nx),  (4) 

and 

R  =  OpHop  +  <rd).  (5) 

Here  kB  is  Boltzmann’s  constant.  T  is  the  lattice  tempera¬ 
ture,  q  is  the  charge  of  a  hole,  c  is  the  dielectric  constant,  is 
the  permittivity  of  free  space,  Nx  and  Nx  are  the  concentra¬ 
tions  of  the  donors  and  traps,  NE  is  the  effective  empty-trap 
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Fig.  1.  Energy-level  diagram  of  band  model  for  photorefraction  in 
BaTiO).  Hole8  are  photoionized  from  level  X  and  recombine  at  X'. 
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Fig.  2.  Crystal  orientation  and  beam  notation  for  two-beam  cou¬ 
pling  in  BaTiO  ,. 


the  relative  contributions  of  photoconductivities  and  dark 
conductivities.  Except  at  low  intensities,  where  the  photo¬ 
conductivity  is  small,  R  is  approximately  equal  to  1.  There¬ 
fore  the  steady-state  photorefractive  effect  is  frequently  de¬ 
scribed  as  being  independent  of  laser  intensity.15-815 
Despite  strong  interest  in  the  applications  of  BaTiO:1  there 
is  a  dearth  of  information  pertaining  to  its  photorefractive 
characteristics  and,  in  particular,  its  intensity-dependent 
effects.  Experiments  are  frequently  conducted  at  intensi¬ 
ties  of  10-100  mW/cm2,  where  the  photoconductivity  domi¬ 
nates  the  dark  conductivity.4  9  Townsend  and  LaMacchia16 
and  later  Rak  et  al.]~  found  an  intensity  dependence  in  the 
photorefractive  effect  in  BaTiO;)  that  saturated  at  approxi¬ 
mately  1  W/cm2,  but  they  did  not  investigate  it  in  any  detail. 
They  attributed  the  intensity  dependence  to  a  diffraction 
efficiency  that  was  proportional  to  the  steady-state  concen¬ 
tration  of  free  carriers.  This  model  is  not  supported  by 
current  photorefractive  theories  and  does  not  explain  other 
aspects  of  the  photorefractive  effect.  More  recently,  Motes 
and  Kim10-11  reported  an  intensity  dependence  in  the  optical 
absorption  coefficient  as  well  as  in  the  two-beam  coupling 
coefficient.  They  showed  that  the  intensity  dependence  in 
the  photorefractive  effect  varied  with  spatial  frequency  and 
was  responsible  for  an  asymmetry  in  the  coupling  coefficient 
with  respect  to  the  c-axis  direction. 


concentration,  and  ad  and  ap  are  the  dark  conductivity  and 
photoconductivity,  respectively.  The  amplitude  of  the 
modulated  refractive  index  n,  is  given  by 

n,  *  (l/2)n3reff£sc,  (6) 

where  n  is  the  index  of  refraction  and  reff  is  the  effective 
electro-optic  coefficient. 

The  grating  formation  in  a  photorefractive  crystal  is  ac¬ 
companied  by  an  energy  redistribution  between  the  two 
interfering  light  beams  that  is  due  to  the  *72  phase  mis¬ 
match  between  the  phase  grating  and  the  light-intensity 
pattern.1213  This  two-beam  coupling  is  shown  schematical¬ 
ly  in  Fig.  2.  The  two  optical  beams  have  transmitted  inten¬ 
sities  Is  and  Ip,  and  the  crystal  c  axis  is  chosen  so  that  the 
signal  beam  (Is)  experiences  gain  defined  as 

Is  (with  IP  on) 

7°*  /s  (with  IP  off)'  (7) 

The  amplitude  gain  can  be  expressed  in  terms  of  exponential 
gain  coefficient12-14  (D: 

(1  +  d0)exp(r  d) 

70  *  i  +  d0exP(rd)  ’  181 

where  d  is  the  interaction  length  and  the  parameter  do  is  the 
intensity  ratio  of  the  two  incident  beams  (do  =  Iso/Ipo)-  The 
exponential  gain  coefficient  is  related  to  the  space-charge 
field  by 

T  =  kn3ref(Esc/m,  (9) 

where  k  is  the  wave  number  of  the  laser  light. 

Equations  tl)-(9)  demonstrate  that  the  exponential  gain 
(two-beam  coupling)  coefficient  is  dependent  on  the  trap 
concentration  and  the  grating  spatial  frequency  but  has  lit¬ 
tle  dependence  on  the  total  light  intensity.  The  only  inten¬ 
sity  dependence  comes  from  the  parameter  R,  which  reflects 


EXPERIMENT 


Although  the  experimental  data  of  Motes  and  Kim11  clearly 
showed  an  intensity  dependence  in  the  absorption  and  pho¬ 
torefractive  two-beam  coupling,  their  photorefractive  data 
were  scattered,  and  comparison  with  theoretical  results  was 
inconclusive.  We  therefore  repeated  those  measurements 
on  a  p- type  BaTiO  (  crystal  that  was  electrically  poled.  The 
BaTiOs  crystal  was  grown  by  Sanders  Associates  and  was 
used  as  grown.  These  nominally  undoped  crystals  typically 
contain  large  concentrations  of  transition-metal  impurities. 
Klein  and  Schwartz5  recently  completed  an  extensive  study 
of  the  impurities  present  in  BaTiO.  (  (grown  by  Sanders)  and 
found  that  iron  was  the  dominant  impurity,  with  concentra¬ 
tions  of  50-150  parts  in  10fi  (ppm). 

The  intensity-dependent  absorption  in  BaTiOi  is  charac¬ 
terized  by  an  absorption  coefficient  that  increases  with  laser 
intensity.  The  intensity-dependent  portion  of  the  absorp¬ 
tion  coefficient,  a/,  can  be  defined  as  the  difference  between 
the  absorption  coefficient  measured  at  high  intensity,  n«. 
and  that  measured  at  low  intensity,  on,  (a;  =  n«  -  a,,).  We 
measured  a/  by  using  the  experimental  arrangement  shown 
in  Fig.  3.  We  used  the  514-nm  line  of  an  argon-ion  laser. 
The  power  of  a  weak  signal  beam  (Is)  was  monitored  by  a 
photodetector.  For  the  absorption  measurements,  the  path 
difference  between  the  pump  (Ip)  and  the  signal  beam  was 
made  much  greater  than  the  coherence  length  to  eliminate 
mutual  photorefractive  effects.  We  also  used  ordinary  po¬ 
larization  to  reduce  the  effects  of  self-induced  photorefrac¬ 
tive  effects.  In  the  presence  of  Ip,  />  was  reduced  owing  to 
the  intensity-dependent  absorption.  We  determined  a/  by 
using  the  relation 


Is  (with  Ip  on) 
Is  (with  Ip  off) 


exp(  —  n/d). 


(101 


where  d  is  interaction  length. 

The  experimental  results  are  shown  in  Fig.  4.  The  varia- 
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514  nm,  S  POLARIZATION 


<3 


Pig.  3.  Experimental  arrangement  for  intensity-dependent  mea¬ 
surements  in  BaTiO,:  BE,  beam  expander;  BS,  beam  splitter; 
DET,  detector;  M's,  mirrors.  For  absorption  measurements,  the 
path  difference  between  Ip  and  Is  was  made  much  greater  than  the 
coherence  length. 


tion  in  the  absorption  coefficient  begins  to  be  observed  at  cw 
intensities  less  than  0.1  W/cm2  and  saturates  at  approxi¬ 
mately  20  W/cm2.  The  intensity-dependent  absorption  ap¬ 
pears  to  be  a  common  characteristic  of  BaTiO  i,  as  it  has  been 
observed  in  all  eight  crystals  that  we  have  tested. 

The  intensity  dependence  of  the  exponential  gain  coeffi¬ 
cient  was  measured  by  using  the  experimental  arrangement 
shown  in  Fig.  3.  Two  coherent  beams  of  an  argon-ion  laser, 
operating  at  514  nm.were  aligned  so  that  the  direction  of  the 
grating  wave  vector  was  parallel  to  the  c  axis  of  the  BaTiO.i 
crystal.  The  beams  were  polarized  become  ordinary  rays. 
The  pump  beam  was  expanded  by  a  2X  beam  expander  to 
ensure  overlap  between  the  two  laser  beams,  and  the  intensi¬ 
ty  of  the  signal  beam  was  much  less  than  the  pump-beam 
intensity  (do  =  0.05).  In  this  geometry  the  signal  beam  Is 
experiences  an  increase  in  intensity  as  described  by  Eqs.  (7) 
and  (8).  When  intensity-dependent  absorption  is  present, 
the  right-hand  side  of  Eq.  (8)  should  be  multiplied  bvthe, 
factor  exp(-oyd)."1  f fcy. 


A  plot  of  the  gain  coefficient  corrected  for  absorption 
versus  grating  period  for  different  pump-beam  intensities  is 
shown  in  Fig.  5.  An  increase  in  beam  intensity  is  seen  to 
result  in  an  increase  in  the  magnitude  of  the  gain  coefficient 
and  a  shift  in  the  peak  of  the  curve  to  smaller  grating  peri¬ 
ods.  Although  an  intensity  dependence  can  be  attributed  to 
the  effects  of  dark  conductivity,  this  cannot  fully  explain  the 
features  observed  in  Fig.  5.  The  intensity-dependent  char¬ 
acteristics  due  to  dark  conductivity  are  manifested  in  Eq. 
(5).  The  coupling  coefficient  should  increase  with  intensity 
from  a  small  value  and  should  level  off  at  intensities  at  which 
the  photoconductivity  dominates  the  dark  conductivity. 
The  small  change  in  the  magnitude  of  T  between  the  two 
lower  intensity  curves  would  indicate  that  the  dark  conduc¬ 
tivity  is  not  dominating  the  photoconductivity  in  this  inten¬ 
sity  range.  The  large  increase  in  T  in  going  to  yet  higher 
intensities  is  not  characteristic  of  a  dark-conductivity  effect. 
Also,  as  seen  from  the  equations  developed  in  the  previous 
section,  the  intensity  dependence  due  to  dark  conductivity 
should  be  independent  of  grating  period.  The  intensity 
dependence  seen  in  Fig.  5  is  more  consistent  with  a  change  in 
the  trap  concentration.  This  point  will  be  discussed  in  more 
detail  below. 

SECONDARY  PHOTOREFRACTIVE  CENTERS 

The  observed  intensity  dependences  in  the  absorption  and 
photorefractive  effect  are  not  in  agreement  with  the  models 
that  are  currently  used  to  describe  photorefraction  in 


S 

>  0.2 


0.01  0.1  1.0  10  100 
INTENSITY  (W/cm2) 

Pig.  4.  Experimental  values  of  intensity-dependent  absorption  co¬ 
efficient  a\  versus  laser  intensity  in  BaTiO;. 
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Fig.  5.  Measured  two-beam  coupling  coefficient  (T)  corrected  for 
absorption  versus  grating  period  for  different  laser  intensities. 
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CB 


Fig.  6.  Schematic  diagram  of  the  secondary-center  model  for 
BaTiOn.  IV i  and  N\'  are  primary  photorefractive  centers,  and  N2  is 
the  secondary  photorefractive  center.  CB,  conduction  band;  VB, 
valence  band. 

BaTiO:i.  In  this  section  we  show  that  the  inclusion  of  sec¬ 
ondary  photorefractive  centers  in  the  band  model  accounts 
for  the  intensity-dependent  effects.18 

The  physical  mechanisms  responsible  for  the  intensity- 
dependent  absorption  in  BaTi(>3  have  not  previously  been 
identified.  Motes  et  al.9  have  shown  that  this  is  not  a  pho¬ 
torefractive  effect.  It  exists  in  unpoled  crystals,  can  be 
induced  at  one  wavelength  and  observed  at  a  second  wave¬ 
length,  and  has  a  rise  time  of  the  induced  absorption  that  is 
an  order  of  magnitude  faster  than  the  grating-formation 
time.  The  effect  is  observed  at  low  intensities,  ruling  out 
two-photon  absorption.  Thermal  effects  are  not  responsi¬ 
ble.  This  is  shown  by  the  decrease  in  absorption  with  in¬ 
creasing  temperature11  and  the  fast  time  response  of  the 
induced  absorption.19  Apparently  either  an  absorption 
mechanism  is  being  created  or  an  existing  mechanism  is 
being  enhanced.  Most  models  used  to  describe  the  photore¬ 
fractive  effect  in  BaTiO-)  consider  only  the  interactions  with 
primary  photorefractive  centers — deep  traps  among  which 
charge  is  redistributed.  These  models  cannot  account  for 
the  intensity-dependent  absorption  at  intensities  of  0.1-10 
W/cm2.  We  suggest  that  the  intensity-dependent  absorp¬ 
tion  is  due  to  secondary  photorefractive  centers  that  exist  in 
addition  to  the  deep  traps.  These  are  traps  that  are  shal¬ 
lower  than  the  primary  centers  and  strongly  ionized.  At  low 
laser  intensity,  there  is,  therefore,  no  measurable  absorption 
due  to  the  secondary  centers.  At  higher  laser  intensities 
these  traps  begin  to  fill  with  holes  generated  in  photoabsorp¬ 
tion  by  the  deep  traps.  Provided  that  the  thermal  ioniza¬ 
tion  rate  is  not  too  high  and  that  the  absorption  cross  section 
is  larger  than  that  for  the  deep  traps,  an  enhancement  in  the 
absorption  coefficient  will  be  observed.  These  secondary 
centers  can  also  be  expected  to  influence  the  photorefractive 
characteristics  in  the  same  intensity  range. 

A  model  that  incorporates  the  secondary  centers  is  shown 
schematically  in  Fig.  6.  This  is  a  modification  of  the  band 
model.  Ni  and  N\  constitute  a  single  set  of  recombination 
centers  (presumably  Fe’u  and  Fe2+,  respectively).  Holes 
are  photoionized  from  N\  and  recombine  at  N\.  The  photo¬ 
ionization  of  Ni  is  taken  to  result  in  the  formation  of  an  N i, 
and,  similarly,  the  hole  recombination  with  an  AY  forms  an 
AY  In  addition  to  these  primary  photorefractive  centers, 
secondary  photorefractive  centers  are  present  and  are  de¬ 
noted  by  N 2.  The  identity  and  characteristics  of  N2  are  yet 
unknown,  and,  for  simplicity,  we  show  only  one  energy  level. 


The  occupation  levels  for  Nt  and  A2  are  described  by  the 
rate  equations 

dNJdt  =  — (S,  /  +  0  l)Nl  +  7  ,(Nir  -  Ni)Nh,  (ID 

dN2/dt  =  -(Sj/  +  82)N2  +  72(A/2T  -  N2)Nh,  (12) 

where  Si  and  S2  are  the  photon-absorption  cross  sections,  71 
and  72  are  the  recombination  coefficients,  di  and  02  are  the 
thermal  ionization  rates,  I  is  the  laser  intensity,  JV2r  and  N\t 
are  the  total  concentration  of  species  N2  and  AY  and  Nh  is 
the  concentration  of  free  carriers  (holes).  (N2t  ~  N2)  is  the 
concentration  of  secondary  photorefractive  centers  that  are 
unoccupied,  and  Nit  =  AY  +  AY  We  include  a  continuity, 
current,  and  Poisson  equation,  respectively,  for  complete¬ 
ness; 

dNh/dt  =  -dA/,/dt  -  dN2/dt  -  V  -j,  (13) 

n  kuT 

j  =  NhfiE  —  B  VN„,  (14) 

<? 

V  •  E  =  (Nh  +  N,  +  N2  -  NA)q/u0,  (15) 

where  m  is  the  hole  mobility,  NA  =  AYO)  +  A/2(0)  +  N/,( 0), 
and  Ni( 0),  N2{0),  and  Nh( 0)  are  the  dark  concentrations. 
N( 0)  reflects  the  dark  conductivity,  and  its  inclusion  pro¬ 
vides  the  associated  intensity  dependence  expressed  in  Eq. 
(5),  as  discussed  above.  Overall  charge  neutrality  is  main¬ 
tained  by  the  presence  of  other  impurities  (e.g.,  oxygen  va¬ 
cancies). 


Fig.  7.  Calculated  intensity-dependent  absorption  coefficient  ver¬ 
sus  laser  intensity. 


Table  1.  Parameter  Values  Used  in  the  Calculations0 


Parameter 

Value 

s, 

2.5  X  10"19  cm-’ 

s2 

5  X  10-’Fcm- 

7l 

5  X  ]0~s  cm:!  sec'1 

y? 

2.5  X  10“°  cm  ’  sec'  ‘ 

d, 

1.5  X  10'4  sec'1 

02 

30  sec'1 

Nit 

2.6  X  10,R  cm':l 

N,(  0) 

2.57  X  101R  cm'1 

N2T 

1.5  X  1017  cm':i 

N2(  0) 

3.8  X  1013  cm'1 

AMO) 

3  X  10s  cm'3 

"For  the  choice  of  the  values  see  the  Discussion  section. 
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Fig.  8.  Calculated  space-charge  field  (Esc)  versus  grating  period 
for  different  intensities. 


To  calculate  the  intensity  dependence  in  the  absorption 
coefficient,  we  assumed  a  uniform  laser  intensity,  set  the 
electric  field  and  current  equal  to  zero,  and  solved  for 
steady-state  conditions.  This  resulted  in  a  set  of  three 
equations  [(11),  (12),  and  (15)]  and  three  unknowns  (Nu  N2, 
and  Nh).  These  equations  did  not  reduce  to  a  convenient 
analytical  form  and  were  solved  numerically.  The  intensi¬ 
ty-dependent  part  of  the  absorption  coefficient  was  then 
given  by 

a,  =  S,[N,  -  N,(0)]  +  S2[N2  -  N2( 0)].  (16) 

Figure  7  shows  the  calculated  steady-state  values  of  a/  as  a 
function  of  laser  intensity  for  the  parameter  values  listed  in 
Table  1.  The  theoretical  results  exhibit  rise  and  saturation 
characteristics  similar  to  the  experimental  data.  This  in¬ 
crease  in  absorption  reflects  the  photoenhanced  population 
of  secondary  centers.  At  low-intensity  illumination  the  sec¬ 
ondary  centers  are  populated  by  capture  of  free  carriers,  but 
thermal  ionization  dominates,  and  there  is  little  change  in 
the  steady-state  concentration  compared  with  dark  condi¬ 
tions.  As  the  intensity  increases,  free  carriers  are  photogen¬ 
erated  from  the  deep  traps  at  a  higher  rate  and  fill  the 
secondary  traps  at  a  rate  comparable  with  the  thermal  ion¬ 
ization  rate.  The  concentration  of  occupied  secondary  cen¬ 
ters  increases;  therefore  the  absorption  coefficient  increases 
because  of  the  high-absorption  cross  section  of  the  second¬ 
ary  centers.  At  yet  higher  intensities,  the  secondary  centers 
begin  to  fill  up,  and  photoionization  limits  the  occupied -trap 
concentration.  The  intensity  region  over  which  this  transi¬ 
tion  takes  place  depends  on  all  the  parameters;  absorption 
cross  sections,  recombination  coefficients,  and,  most  criti¬ 
cally,  thermal  ionization  rate. 

To  calculate  the  steady-state  space-charge  field  we  used  a 
modulated  intensity  distribution  given  by  Eq.  (1).  We  as¬ 
sumed  that  m  «  1  and  solved  for  the  lowest  Fourier  compo¬ 
nent  of  the  grating  by  writing  the  variables  as  zeroth-  and 
first-order  terms: 


N ,  =  N10  +  Nu  exp (iKz), 

(17) 

N2  =  N20  +  N21  exp  (iKz), 

(18) 

-  Nh0  +  Nhl  exp(iKz), 

(19) 

E  =  E0  +  Esc  exp(iKz), 

(20) 

)  ~  Jo  +  h  exp  (iKz). 

(21) 

Eo  and  Esc  correspond  to  the  applied  electric  field  and  in¬ 
duced  space-charge  field,  respectively.  We  considered  only 
the  diffusion  problem  where  Eo  and  jo  were  zero.  Solution 
for  the  steady  state  proceeded  by  substituting  Eqs.  (17)— (21) 
into  Eqs.  (11)— (15)  and  keeping  only  zeroth-  and  first-order 
terms.  This  yielded  a  set  of  coupled  nonlinear  equations  in 
the  zeroth-order  terms  that  were  solved  numerically  for  N 10, 
N20,  and  Nho,  which  were  then  used  as  parameters  for  the  set 
of  linear  equations  in  the  first-order  terms.  These  equa¬ 
tions  were  then  solved  numerically  for  Esc  as  well  as  for  Nu, 
N2i,  and  Nfci. 

Figure  8  shows  the  calculated  dependence  of  Esc  on  grat¬ 
ing  period  and  intensity  for  the  parameter  values  listed  in 
Table  1.  At  low  laser  intensity  few  secondary  centers  are 
occupied,  and  the  secondary-center  model  gives  the  same 
result  as  the  Kukhtarev  model.  As  the  intensity  increases, 
the  steady-state  concentration  of  secondary  centers  (N 2) 
also  increases.  This  in  turn  results  in  an  increase  in  the 
effective  empty-trap  concentration.  Consequently  the 
magnitude  of  the  space-charge  field  is  enhanced,  and  the 
peak  is  shifted  to  smaller  grating  periods.  Although  the 
secondary  centers  are  traps,  which  are  present  even  at  low 
intensities,  they  do  not  contribute  to  the  space-charge  field 
until  they  are  occupied  with  charge.  Trapping  of  a  hole  by  a 
secondary  center  also  influences  the  concentration  distribu¬ 
tion  of  primary  centers.  If  secondary  centers  are  not 
present,  the  steady-state  concentrations  of  Ni  and  Nj  re¬ 
main  essentially  unchanged  from  their  initial  values.  But 
each  hole  trapped  by  a  secondary  center  results  in  a  decrease 
in  the  Ni  concentration  and  an  increase  in  the  A)] 'concentra¬ 
tion.  Therefore  the  increase  in  the  effective  empty-trap 
concentration  is  given  approximately  by  2[A^2  -  N2(0)]. 

The  effect  of  the  secondary  centers  is  more  complicated 
than  a  simple  increase  in  the  number  of  participating  traps. 
For  a  given  grating  period  the  magnitude  of  Esc  is  propor¬ 
tional  to  the  magnitude  of  the  modulation  in  the  ionic  space 
charge  (i.e.,  Nn  +  N21).  Although  the  modulation  of  N\  was 
always  dephased  180°  with  respect  to  the  light-intensity 
sinusoid,  the  modulation  in  N 2  was  found  to  be  either  in 
phase  or  180°  out  of  phase,  depending  on  intensity  and 
grating  period.  In  general,  for  low  intensities  N2  was  in 
phase  and  at  higher  intensities  and  shorter  grating  periods 
N2  was  out  of  phase  with  the  light  intensity.  For  the  results 
shown  in  Fig.  8,  N2  was  in  phase  for  the  two  lower-intensity 
curves,  and  at  2  W/cm2,  N2  was  out  of  phase  for  grating 
periods  of  less  than  0.36  ftm.  The  enhanced  space-charge 
field  at  small  grating  periods  was  due  primarily  to  an  in¬ 
crease  in  Nn,  which  dominated  N2].  As  the  grating  period 
increased,  the  magnitude  of  Nu  decreased,  whereas  N2i  in¬ 
creased.  At  a  grating  period  of  1.65  nm,  the  space-charge 
field  was  only  slightly  enhanced  at  an  intensity  of  2  W/cm2 
over  that  at  0.002  W/cm2,  despite  an  increase  in  the  value  of 
Nu  by  a  factor  of  4.  This  was  because  the  increase  in  Nn 
was  offset  by  the  modulation  in  N2.  We  applied  the  model 
to  other  parameter  values  and  intensities  but  found  no  com¬ 
bination  that  allowed  N2i  to  become  greater  than  Nu  and 
thereby  to  reverse  the  sign  of  the  coupling. 

DISCUSSION 

Values  for  the  various  parameters  used  in  this  model  are  not 
well  known,  particularly  for  the  secondary  centers.  The 
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values  listed  in  Table  1  were  chosen  to  be  consistent  with  our 
data  and  the  current  understanding  of  BaTi03,  as  discussed 
below.  These  values  do  not  necessarily  represent  a  best  fit 
to  the  data.  There  are  too  many  parameters  for  which  more- 
detailed  knowledge  is  required  to  justify  such  an  analysis. 

Iron  impurities  are  thought  to  be  the  primary  photorefrac- 
tive  centers  in  BaTi03  and  are  usually  present  in  concentra¬ 
tions  of  approximately  100  ppm  (5  X  1018  cm'3),  with  an 
Fe2+  to-Fe3+  ratio  of  approximately  0.01  for  p- type  crys¬ 
tals.8  We  estimated  the  N ('  (Fe2+)  concentration  by  using 
the  effective  empty-trap  concentration  of  2.6  X  1016  cm-3,  as 
determined  from  the  two-beam  coupling  data  at  low  intensi¬ 
ty.  We  then  took  N i(0)  to  be  100  times  this  value.  Klein 
and  Schwartz5  determined  impurity  concentrations  and  ab¬ 
sorption  coefficients  at  442  nm  for  several  BaTi03  crystals. 
From  their  data  we  estimated  the  absorption  cross  section 
for  N\  to  be  approximately  2.5  X  10-19  cm2  at  514  nm.  The 
recombination  coefficient  is  thought  to  be  approximately  5 
X  10-8  cm3/sec,  consistent  with  a  free-carrier  lifetime  of 
approximately  1  nsec.8  The  deep  traps  have  a  low  thermal 
ionization  rate,  which  is  reflected  in  long  grating-storage 
times.  We  used  the  value  1.5  X  10-4  sec-1  to  be  consistent 
with  a  finite  dark  level  of  carriers. 

Characteristics  of  the  secondary  centers  can  be  obtained 
through  interpretation  of  the  intensity-dependent  absorp¬ 
tion  and  photorefractive  data.  The  increase  in  the  effective 
empty-trap  density  inferred  from  the  two-beam  coupling 
data  implies  an  increase  in  the  occupied  secondary-center 
concentration  of  approximately  4  X  1016  cm-3  at  2  W/cm2. 
This  can  be  combined  with  the  intensity-dependent  absorp¬ 
tion  data  to  give  an  absorption  cross  section  of  approximate¬ 
ly  5.25  X  10-18  cm2.  The  increase  in  occupied-secondary- 
center  concentration  saturated  at  approximately  6  X  1016 
cm'9.  This  suggests  that  the  total  concentration  of  second¬ 
ary  centers  (N2t)  was  approximately  1  X  1017  to  2  X  1017 
cm-3.  We  have  estimated  the  thermal-ionization  rate  to  be 
between  20  and  100  sec-1,  based  on  the  temporal  decay  of  the 
intensity-dependent  absoprtion.  This  provides  only  an  ap¬ 
proximate  value  because  the  measured  decay  time  is  strong¬ 
ly  influenced  by  the  recombination  coefficients  and  concen¬ 
trations  of  N\  and  N2.  We  find  that  a  value  between  20  and 
40  sec-1  best  reflects  our  absorption  data. 

The  secondary  center,  N2,  is  thus  characterized  as  a  highly 
ionized  trap  with  a  moderate  thermal -ionization  rate.  The 
presence  of  such  an  impurity  is  not  in  conflict  with  the 
known  defect  chemistry  of  BaTi03.  As-grown  crystals  are 
generally  p  type  and  are  characterized  by  a  high  dark  resis¬ 
tivity  of  10u-1016  f!  cm  and  therefore  a  low  free-carrier 
concentration  in  the  dark.20-22  Although  BaTi03  normally 
contains  a  high  concentration  of  transition-metal  impurities 
and  barium  vacancies  that  act  as  acceptors,  a  high  concen¬ 
tration  of  oxygen  vacancies  is  usually  present,  which  com¬ 
pensates  for  the  acceptor  level,  even  when  grown  in  a  slightly 
oxidizing  environment.5-20 

We  have  not  included  photoabsorption  by  the  AV  impuri¬ 
ty  states.  Besides  providing  centers  for  hole  recombination, 
Ni'  may  also  provide  a  center  for  electron  photoexcitation. 
This  process  has  recently  been  considered  as  a  source  for 
electron-hole  competition  in  photorefraction.23-24  In  p-type 
BaTi03,  in  the  absence  of  secondary  traps,  this  absorption 
mechanism  has  negligible  effect  on  the  intensity -dependent 
absorption  at  the  low  laser  intensities  considered  here. 


Much-higher  intensities  would  be  required  to  change  the 
concentrations  of  N\  and  N However,  as  discussed  above, 
when  secondary  centers  are  present,  the  concentration  of  Nj 
is  increased  along  with  the  concentration  of  occupied  sec¬ 
ondary  traps.  By  neglecting  absorption  due  to  N T,  we  have 
effectively  assumed  that  the  absorption  cross  section  Si  is 
small  compared  with  S2  and  have  attributed  all  the  en¬ 
hanced  absorption  to  the  secondary  centers.  If  absorption 
by  Ni  is  important,  then  our  conclusions  about  the  charac¬ 
teristics  of  the  secondary  centers,  such  as  the  value  of  the 
recombination  coefficient  and  the  absorption  cross  section, 
would  need  to  be  modified  accordingly.  However,  the  es¬ 
sential  feature  of  the  model  would  be  unaffected.  The  num¬ 
ber  of  secondary  centers  involved,  as  indicated  by  the  pho¬ 
torefractive  data,  would  still  be  the  same.  Secondary  traps 
would  still  be  necessary  to  contribute  to  charge  redistribu¬ 
tion  in  either  case.  We  have  examined  the  case  of  electron 
traps  populated  through  absorption  by  Ni Although  these 
traps  could  account  for  the  intensity-dependent  absorption, 
the  secondary  hole  traps  were  still  necessary  to  account  for 
the  intensity-dependent  photorefractive  effects. 

The  secondary-center  model  presented  here  was  intended 
to  reflect  conditions  appropriate  for  p- type  BaTi03.  As 
such,  the  calculations  were  limited  to  hole  production  and  a 
single-carrier  photorefractive  model.  Our  objective  was  to 
provide  a  basis  for  understanding  the  intensity-dependent 
effects  in  BaTi03.  Extension  of  the  model  to  include  multi¬ 
ple  secondary  traps,  both  electron  and  hole,  and  mixed  con¬ 
ductivity  is  straightforward.  However,  extension  to  a  more 
elaborate  model  would  also  require  additional  parameters 
for  which  there  is  little  or  no  knowledge. 

CONCLUSION 

We  have  attributed  the  intensity-dependent  absorption  and 
photorefractive  effects  in  BaTi03  at  cw  laser  intensities  to 
the  presence  of  secondary  photorefractive  centers.  These 
are  intermediate-level  traps  that  are  not  completely  occu¬ 
pied  at  room  temperature  and  compete  for  charge  carriers 
photoionized  from  the  deep  traps.  The  increase  in  the  occu¬ 
pied  concentration  with  higher  laser  intensity  results  in  an 
enhanced  absorption  coefficient.  The  secondary  centers 
also  play  a  role  in  photorefraction  by  providing  an  intensity- 
dependent  trap  concentration. 
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